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ABSTRACT 



From a study of low-pass filtered hourly sea levels for Sattahip 
and Ko-Lak in the northern Gulf of Thailand during the period 1960 
through 1966, the sea level variations due to meteorological effects 
are found. The mean annual variation of the filtered sea level, 
which averages 0.5 meters, is consistent with the climatological mean 
annual wind variations. This sea level change appears to be a 
response to the Ekman transport which would be expected from the 
seasonal monsoon wind regimes. Sea level response to atmospheric 
pressure is negligible compared with response to wind. Analyses 
performed to find relations between the filtered sea level and the 
gulf-wide geostrophic wind and local surface wind show that the sea 
level slopes upward across the gulf in the same direction as the 
local wind blows, the response being coherent in the frequency band 
0.083-0.117 cycles per day (period 8 to 12 days). This wind set- 
up effect is clearly secondary to the Ekman transport in inducing 
the seasonal sea level variations observed in the gulf. 
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I. INTRODUCTION 



In 1934 the Hydrographic Department of the Royal Thai Navy installed 
a tide gauge to record sea level at the mouth of the Chao-Phaya river 
in the Gulf of Thailand and sent the data to the U.S. Coast and 
Geodetic Survey (U.S. A.) for analysis and prediction. In 1959, 
the Royal Thai Navy bought the Doodson-Lege' prediction machine from 
England and used it to predict tides in 1960. It also set up tide 
gauges at many stations along the coast of the gulf until today it 
is possible to predict tides at about 29 stations. Since 1972, tide 
prediction has been done with the computing machine (I.B.M.360) at the 
Supreme Command, Ministry of Defence, but still the Doodson tide prediction 
machine is used also for comparison. For the analysis by the Doodson 
harmonic method (Doodson, 1941) about 30 constituents are calculated. 

Tide in the Gulf of Thailand is the shallow tide; the coast and 
continental shelf are irregular. This gulf is bordered by Thailand on 

the north, Cambodia and Vietnam on the east, and Malaysia on the west, 
and its mouth opens to the South China Sea (Figure 1). It has a width 
of about 200 nm., length about 400 nm., and average depth about 80 meters. 
At the northern end, the head of the gulf narrows. This region is about 
60 nm. wide and about 80 nm. long, with average depth about 20-30 meters. 
The gulf is affected by monsoon winds and by tropical storms. There are 
3 seasons: 1) southwest monsoon, April -September (about 6 months); 2) 

northeast monsoon, October-January (about 4 months); 3) southeast season, 
February-March (about 2 months) with the wind from the South China Sea 
which has high humidity. 
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The analysis and prediction done by the Royal Thai Navy treats only 
the forcing with periods of the astronomical tides. This thesis will 
be concerned with seasonal variation and other largely non-astronomical 
effects on sea level, and their causal relationship to meteorological 
conditions. It will be shown that non-astronomical events contribute 
substantially to the water level variations in the gulf and are therefore 
important in prediction of sea level. 
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II. NON-ASTRONOMICAL FLUCTUATION IN SEA LEVEL 



The observed elevation of the ocean's surface at any location and 
time is influenced by a combination of factors, including the large scale 
atmospheric circulation. Generally, the 12.4- and 24.8-hour tide-raising 
forces of the earth-moon-sun system generate the most conspicuous and 
easily recognized of all ocean tides. The observed sea level, however, 
is modified by combinations of other less predictable and quantifiable 
factors, operating individually or in concert. "These include atmos- 
pherically driven on-and offshore Ekman mass transport, direct set-up 
or set-down of water due to winds, spatial and temporal variations in 
atmospheric pressure, wind generated waves and swell, density variations 
in the water column, crustal subsidence or uplift, eustatic sea level 
variations resulting from changes in the mass of the polar ice packs, 
surge from distant storms or tsunamis, and low frequency astronomical 
tide-raising forces". (Chelton, 1980). 

In 1969, the Fleet Numerical Weather Central (now called "Fleet 
Numerical Oceanography Center" or "FNOC"), Monterey, California (U.S.A.) 
studied tides in the Gulf of Thailand, and pointed out some of the features 
of the atmospheric overtides. According to FNOC (1969): 

"a) Any wind can cause changes of sea level in addition to the 
astronomical tides. Thus, there is a difference in mean sea level during 
different monsoon winds, the difference varying from location to location 
and being greatly influenced by the bathymetry at a given location. 
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b) There is a considerable mass transport of water by waves in 
shallow water. This mass transport can cause strong alongshore currents, 
which might also, during storms, cause a considerable shift of the sand 
bank. If there are sand bars or coral reefs off the coast, a considerable 
amount of water is carried over them by waves. This wave "set-up" can 
account, occasionally, for as much as three-foot changes of water level. 

c) The sea level also acts as an inverted barometer. A 1-mb 
decrease in atmospheric pressure will cause a 1-cm. rise in the sea level. 

d) If the atmospheric pressure disturbance moves with a speed 
comparable to the speed of a shallow water wave, v =7gT; g = acceleration 
of gravity, and d - depth in meters; the atmospheric pressure effects 
are overamplified. 

e) Twelve inches or more of rainfall during a 24-hour period can 

be caused by typhoons. This river run-off can cause water level changes 
in addition to those due to wind and atmospheric factors. 

f) There are relatively great differences in the storm surge chara- 
cteristics between geographical points only 10 miles or so apart. 

g) There are "overshooting" or "resurge" effects, especially in 
bays and estuaries in connection with storm surges. 

h) In extensive shallow water areas the local wind effects on the 
water and movement of the water within the bay are as important as the 
external surges entering the area". 

It is obvious that improved sea level prediction can be the result 
of applying the non-astronomical variation to that given by the astrono- 
mical tides alone. One purpose of thesis is to assess the importance 
of some of the non-astronomical factors in the Gulf of Thailand. 



14 



III. DATA 



Sea level data and mean annual climatological elements from Sattahip 
and Ko-Lak, and daily atmospheric pressure and wind from Bangkok, Sattahip, 
Singapore, and Saigon were used in this study. 

The hourly tide heights from Sattahip and Ko-Lak were provided by 
the Hydrographic Department of the Royal Thai Navy. Sattahip station 
is on the east coast and Ko-Lak station is on the west coast in the northern 
part of the gulf. They are 80 nm. apart and situated as shown in Figure 
1. Station Sattahip is in Sattahip Bay, which has dimensions 4 nm. 
in width, 1 nm. in length, and average depth about 6 meters. 

About 20 years of hourly readings were available for both stations 
and these have been tabulated in the standard table form (Form 362, 

U.S. Coast and Geodetic Survey) for the period 1 January 1960 through 
31 December 1979. Some of these data were punched on computer cards 
for use in later analyses. Because it takes a lot of time to punch the 
cards and because a shorter time is enough for the present analysis, 
only 7 years of data, for the period 1 January 1960 through 31 December 
1966, were used. The hourly heights had numerous short gaps and some 
long gaps (Table 1), the longest being 14 days in 1960 at Sattahip station. 
These gaps were filled by conventional harmonic prediction. 

The mean monthly climatological data presented herein from Sattahip 
and Ko-Lak were monthly averages from the years 1951 through 1975. 

These data were provided by the Meteorological Department, Ministry of 
communications, Bangkok, Thailand. 
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Daily atmospheric pressure and surface wind data were obtained from 
the NAVOCEANCOMDET, Asheville, N.C. These data were available for 
Bangkok, Sattahip, Singapore, and Saigon during January 1 960-December 
1965. The data were recorded usually for the hours 00, 03, 06, 09, 12, 
and 18 GCT. These records, however, are not regularly available for the 
same hour for every day, except for the time 0600 hours in 1962 which 
has continuous values. Sattahip has a lot of missing data, so it was 
not considered in this analysis. Accordingly, meteorological data at 
0600 hours in 1962 for Bangkok, Saigon, and Singapore were used. 
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TABLE 1 



Gaps in Sea Level Records 



Station Ko-Lak 



0000 


8 Mar. 


to 


1100 


10 Mar. 


1960 


(2 days) 


1600 


26 Oct. 


to 


0600 


29 Oct. 


1960 


(4 days) 


2000 


10 Dec. 


to 


0600 


13 Dec. 


1960 


(3 days) 


1200 


6 Jul . 


to 


1000 


14 Jul . 


1961 


(9 days) 


0000 


15 Jul. 


to 


1500 


16 Jul . 


1962 


(2 days) 


0900 


16 Jun . 


to 


1500 


19 Jun. 


1963 


(4 days) 


1800 


29 Aug. 


to 


2300 


31 Aug. 


1963 


(2 days) 


1600 


12 Dec. 


to 


2300 


21 Dec. 


1964 


(9 days) 


0900 


25 Nov. 


to 


2300 


30 Nov. 


1965 


(5 days) 


0000 


1 Dec. 


to 


0700 


6 Dec. 


1965 


(5 days) 


0000 


1 Nov . 


to 


1600 


4 Nov. 


1966 


(3 days) 


0000 


1 Dec . 


to 


0800 


1 Dec. 


1966 


(8 hours) 


Station Sattahip 












2000 


5 Sep. 


to 


1600 


20 Sep. 


1960 


(14 days) 


0400 


16 Nov. 


to 


0800 


23 Nov. 


1960 


(8 days) 


1800 


6 Dec. 


to 


1600 


10 Dec. 


1960 


(4 days) 


1100 


31 May 


to 


2300 


31 May 


1965 


(13 hours) 


0000 


20 Mar. 


to 


1000 


22 Mar. 


1966 


(3 days) 


1100 


24 Apr. 


to 


1000 


27 Apr. 


1966 


(3 days) 


0000 


6 Jun . 


to 


0900 


8 Jun . 


1966 


(3 days) 


1800 


20 Jun . 


to 


2300 


20 Jun . 


1966 


(6 hours) 


1400 


30 Oct. 


to 


2300 


31 Oct. 


1966 


(2 days) 


0000 


1 Nov. 


to 


1000 


2 Nov . 


1966 


(2 days) 


1100 


1 1 Dec . 


to 


1000 


14 Dec. 


1966 


(3 days) 
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IV. ANALYSIS OF SEA LEVEL HEIGHTS 



The hourly heights for both stations, Ko-Lak and Sattahip, can be 
reduced to approximately the level of the non-astronomical water level 
by separating the prominent semi-diurnal and diurnal tide components from 
the non-harmonic components of the hourly height. Then the non-astrono- 
mical part of the height change can be seen clearly. This process, a 
low-pass filtering, was done by computing running means of the hourly 
heights using averaging intervals determined by the period of harmonic 
components to be suppressed. Godin's (1966) filter accomplished this 
effectively. The filter performs an unweighted arithmetic running mean 
of the data set, passing over the data three times. It averages sequences 
of 24 values on the first and second passes, and 25 values on the third- 
pass. Appendix A shows the unfiltered hourly heights, the filtered semi- 
diurnal and diurnal tides (eta(o)-eta( f) ) , and the filtered residual 
hourly heights, with each figure containing one year of data. 

The filtered heights from both tide stations show a lot of noise 
throughout the year, much of which is at clearly non-astronomical frequ- 
encies. The maximum difference between high and low water is 8-10 dm. 
in all the years. Large and rapid changes at each station occur at almost 
the same times and are of similar height. The spikes (both maxima and 
minima) are about 5-6 dm. in height over intervals of 12-30 days. 

These spikes usually are largest during the period January through March. 
These kinds of variations might be expected as the result of meteoro- 
logical disturbances, which have similar time scales. 
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To study the sea level and its relation to meteorological forcing 
more closely, spectrum and cross-spectrum analyses of data were made. 

Two problems that arise in spectrum analysis are aliasing from high 
frequencies and leakage, usually from low frequency spectral peaks. 

These problems are discussed below. 

A. ALIASING 

If energy is present at frequencies higher than the highest analy- 
zable (or Nyquist) frequency (1/2 A), this energy will appear at an apparent 
frequency lower than the actual frequency, where is the sampling interval. 
This energy is said to be "aliased" into the lower frequency which is 
also lower than 1/2A. 

The phenomenon of aliasing is important in practical data analysis. 

The sampling interval (^) must be short enough to cover the full frequency 
range of the continuous time series. Otherwise the spectrum derived 
from equally spaced samples will differ from the true spectrum because 
of aliasing. In some cases the only way to be certain that this condi- 
tion is met may be to filter the time series to remove intentionally all 
frequency components higher than 1/2 a before beginning the analysis 
(.Newland, 1974). This procedure is of course impossible if the data 
set is available only at the sampling interval (A). 

One source of aliasing in sea level data in the present analysis 
may be a resonant frequency in the Gulf of Thailand, a frequency which 
is estimated to be close to the Nyquist frequency (f = 0.5 cpd) for the 
once daily sea-level heights left after the use of Godin's filter. This 
frequency can be estimated as follows. 
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Since the Gulf of Thailand is shallow, the resonant period can be 



estimated from the Sverdrup wave phase speed. 




where 

g = gravity acceleration 
h = the mean depth 
6 = radian frequency 
f = coriolis parameter 

Resonant fundamental frequency for an open rectangular basin of 
constant depth, fr, is 
fr = c/41 

where 1 is the length of the gulf. Assuming the gulf is rectangular 
with a mean depth (h) of 60 meters, length (1) of 400 nm., and mean 
latitude of 10 degrees, the estimated resonant frequency is 0.6 cycles 
per day (cpd) . 

For most spectra and cross-spectra involving the meteorological data 
used here, the sampling interval is A = 1 day. So the Nyquist frequency 
(fn) = 1/2A = 0.5 cpd. This means that, if the estimated free 
(resonant) frequency is correct, 
fn < fr 

Thus an important frequency may be aliased into the frequency 
f = 2fn - fr. 

In fact, the estimate of resonant frequency is crude. So one can 
say that the resonant frequency is either barely resolved or not quite 
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resolved. As a result, energy at the resonant frequency may account 
for spectral peaks seen near the Nyquist frequency or at the Nyquist 
frequency itself. 

B. LEAKAGE 

Another problem in spectrum analysis is "leakage" from a low frequency 
or other spectral peak. This can be suppressed by prewhitening. 

According to Newland (1974), "The process of 'prewhitening' involves 
operating on raw data before spectral analysis in such a way that the 
processed data has an approximately flat spectrum. Normalizing to 
zero mean and correcting for slow trends are two examples of this. 

Such initial processing of data before analysis is carried out for two 
reasons. One is that, because large peaks in the spectrum are reduced, 
spectral 'leakage' due to the imperfect shape of the spectral window is 
reduced. The other arises because the standard error of a measurement 
is proportional to the mean value of this measurement. The variance 
of a large spectral estimate is therefore proportionally more than the 
variance of a small estimate. By whitening the spectrum before analysis, 
the estimated spectrum can be given uniform accuracy at all frequencies". 

In order to test the spectrum analysis procedure as to the need for 
prewhitening of the low-pass filtered sea levels at Sattahip and Ko- 
Lak, spectra from the 1962 data were calculated both with and without 
prewhitening by the program FTFREQ (IMSL, 1980). These spectra are 
shown in Figure 2. The spectra with prewhitening were obtained by 
replacing the ith term of each series Xi was replaced by Xi - Xi-1 
before analysis. The two processes (with and without prewhitening) 
covered the 6-month periods (October-March and April -September) . 
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The maximum number of lags in cross-and autocorrelation was m = 30, the 
number of terms was n = 180, and the number of degrees of freedom for the 
estimation of confidence limit was d.f. = 2n/m = 12. 

The spectra from both analyses are very nearly the same. This means 
that leakage is not a problem and that one can analyze similar data 
without prewhitening. 
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Figure 2. Spectra of Filtered Sea Level (1962) 
(a) with prewhitening (b) without prewhitening 
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V. 



RESULTS OF THE ANALYSIS 



A. SEASONAL VARIATION OF SEA LEVEL 

To examine the seasonal variation free of the short period variation, 
additional smoothing was done of the residual hourly sea levels obtained 
from application of Godin's filter. A 28-day period (672 hours) was 
selected because it minimizes the prominent fortnightly cycle and the 
bi -fortnightly variation of spring/neap tides, it suppress the relatively 
high-frequency non-astronomical noise, it is the length of a lunar month, 
and it approximates the commonly used comparison period of one month. 

The time series of Appendix B show the hourly heights treated by 
Godin's filter (24-24-25 hours) and the 28-day filter (672 hours). 

The figures show that sea level usually increases from July to January 
and decreases from February to June at both stations. This corresponds 
to increasing sea level during the northeast monsoon, and decreasing sea 
level during the southwest monsoon. In these figures the sea level 
still shows great variability. These long-period fluctuations may be 
caused by meteorological and oceanographic effects and perhaps also by 
local bathymetric changes, such as formation and removal of sand bars. 

B. MEAN ANNUAL VARIATION OF SEA LEVEL 

The relation between seasonal sea level and atmospheric changes can 
be studied by comparing the annual sea level change with annual atmospheric 
pressure changes and variations in the local wind. The mean annual 
variation of meteorological elements is given in Table 2. The mean 
monthly sea levels were computed by averaging monthly levels from 1960 
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to 1966. The monthly pressure and wind data are averages from 1951 

to 1975. These values are plotted in Figure 3. It is seen that sea 

level decreases from February to June and increases from June to February, 
which is in general agreement with the seasonal tendency of the mean 
monthly atmospheric pressure. The effect of atmospheric pressure on 
sea level can be estimated by applying the inverted barometer (hydrostatic) 
relationship, i.e., sea level is depressed by approximately 1 centimeter 
for every millibar of pressure increase. The mean annual variation of 
sea level, however, is seen in Figure 3 to be about 8 times the inverse 
barometer effect, and in addition the observed relationship is direct 
rather than inverse. This shows that the annual variation of atmospheric 
pressure accounts for very little of the annual sea level variation in 
this region. 

Now consider the relations between the annual wind variation and the 
annual variation of sea levels at Sattahip and Ko-Lak, these variations 
being shown in Table 2. and graphically in Figure 3. These variations 
can be viewed as having two components, one representing the mean rise 
in level at the northern end of the gulf, the other representing the 
cross-gulf slope between the Ko-Lak and Sattahip. It is obvious from 
the similarity between the annual regimes at the two stations (Figure 3c) 
that the variations in mean level are much greater than those associated 
with cross-gulf slopes. The variation in mean sea level is consistent, 
qualitatively, with a tendency for water to leave the gulf by Ekman 
transport associated with winds from the south (as in February through 
June) and to return to the gulf when the winds are from the north (as 
in October through January). The cross-shelf slopes are small, but there 
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is a tendency for Sattahip to be relatively high during the period of west 
or southwest winds (June through September). This is consistent with 
a shallow water set-up response in phase with those winds. 

C. SPECTRUM OF SEA LEVEL 

To examine the sea level and its relation to forcing processes more 
closely, spectrum analysis of the data was used. The spectrum of sea 
level represents the distribution of "energy", or mean-square oscillation, 
over a range of frequencies. The frequency components present in a 
particular time series can be found and their significance evaluated. 

The method can also be applied to studying the relation between two time 
series, for example, the sea level at two neighboring stations for the 
same period of time. According to Hamon (1962), "In such cases the 
degree of association ('coherence') is obtained as a function of frequency; 
and the time displacement of one series relative to the other is found 
in terms of the 'phase difference', also as a function of frequency. 

The coherence and phase difference can give useful information even when 
the spectrums of the separate time series are relatively featureless". 

Characteristic features of the sea level spectrum can be considered 
in three frequency regions according to Godin (1966): 

"1) The low frequency band, defined over an interval extending from 
0 to 1 cycle per day; 

2) The central band, limited to 1 to 6 cycles per day; 

3) Finally, the high frequency band containing all the frequencies 
above 6 cycles per day". 
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Table 2. Mean monthly sea level, pressure, and wind 
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The central band contains the diurnal and semi-diurnal periods and 
is the most energetic part of the spectrum. The low-frequency band 

contains the seasonal cycle, yearly variations, and meteorological 
disturbances. 

In this analysis, the variance densities of the hourly heights at 
Sattahip and Ko-Lak were calculated, and the coherence, phase, and cross- 
spectra of the two data sets were obtained. These spectra are examined 
in two frequency regions; first, including the central band (0.04 to 
0.25 cycles per hour) and second, in the low frequency band (0 to 0.04 
cycles per hour) with the central band eliminated by the low pass filter 
(Godin, 1966). Both analyses were performed on three-month blocks of 
hourly height ( January-March, April-June, etc.,) for seven years of data. 
These analyses are graphically presented in Appendix C and Appendix D. 

The spectrum for each data set has resolution f = 1/512 hr”^ and 8 
degrees of freedom. 

The spectra of the raw hourly heights, shown in Appendix C, which are 
dominated by frequencies lying in the central band, show the prominent 
peaks of the diurnal and semidiurnal tidal components at frequencies about 
0.041 and 0.08 cycles per hour. At these frequencies, the coherence 
between the two stations tends toward unity, which means that they are 
fully related. The phase graph shows that these tidal components at 
Ko-Lak lag behind those at Sattahip by 3-20 degrees, which corresponds 
0.2-1 .7 hours. The amplitude of the diurnal tide is much higher than 
the semi-diurnal tide, which indicates that usually the tide has diurnal 
form. 
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To consider the seasonal variation and others of low frequency, it 
is useful to eliminate the effect of the prominant semi-diurnal and diurnal 
astronomical tides. The spectra are made free of these astronomical tides 
by low pass filtering (Godin, 1966). This essentially filters out high- 
frequency fluctuations. In addition, the linear trend was removed from 
the records to make the time series as stationary as possible (Roden, 1966). 
The results are graphically presented in Appendix D and summarized for 
each three-month period for seven years in Figures 4 and 5. These 
spectra show a large interannual (year-to-year) variation, especially in 
the January-March period. The years 1962 and 1964 have 8 or 10 times the 

energy of 1960 and 1965. This means that the non-astronomical sea level 

variations can not be predicted well from a mean annual variation. It 
is important to know the variation from year-to-year. 

To increase stability of the spectral estimates and to show the mean 
spectral form, the spectra for each three-month period in Figures 4 and 
5, were ensemble averaged over the seven-year period. These averaged 
spectra are presented in Figures 6 and 7; they have resolution f = 

1/512 hr~^ and 56 degrees of freedom. The energy in frequencies above 
0.021 cycles per hour quickly tends toward zero, and consequently these 
frequencies are omitted. This cut-off held true for both raw and filtered 
hourly data, with the exception of the extremely strong energy peaks 
at the 12.4 and 24.8 hr. lunar-solar tidal periods in the hourly data. 

The spectra show the highest energy in the period January-March at both 
stations. This shows that the large variations in the sea level occur 
in the period of the northeast monsoon season. 
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To study the interannual variation for the entire year, the spectra 
for each three-month block are ensemble averaged over each whole year 
(as in Figures 6 and 7; resolution f = 1/512 hr "* and 56 degrees of 
freedom). These spectra are shown in Figure 8; they have the same resolution 
and 224 degrees of freedom. The spectra show high peaks at frequency 
0.00195 cycles per hour (period 21.4 days) and rapid decrease to frequency 
0.0058 cycles per hour (period 7.1 days). For these spectra the energy 
at frequency 1/512 hr~^ has contributions from the annual frequency and 
other low frequencies. Their energy all appears at that (first resolvable) 
frequency. This is the characteristic feature of a "red spectrum", 
which defines the increase in energy density with decreasing frequency. 

The spectra of non-periodic sea level fluctuations are decidedly red, 
showing a concentration of spectral energy at low frequencies (Roden, 1960). 

In Figure 8, the spectra are seen to be very similar at the two stations 
(Ko-Lak and Sattahip). This suggests that the low-frequency non-astro- 
nomical components of sea level at both stations may be responding to 
the same driving mechanisms, perhaps variations in the monsoon wind. 

For the analysis of seasonal differences the spectra from each station 
were ensemble averaged over each monsoon season, and are shown in Figure 
9. The spectra in each three-month period in Figures 6 and 7, (with 
resolution f = 1/51 2 - ^ and 56 degrees of freedom) were then ensemble 
averaged in 6-month periods in order to include the northeast monsoon 
(October-March) and southwest monsoon (April -September) seasons. The 
resulting spectra have the same resolution and 112 degrees of freedom. 

Due to the constraints imposed by the three-month analysis process, the 
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Figure 4. Seven years of filtered sea-level spectra for three- 
month periods at Sattahip C8 degrees of freedom). 
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Figure 5. Seven years of 
month periods at Ko-Lak (8 



filtered sea-level spectra 
degrees of freedom) . 



for three- 
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Figure 6. Seven-year ensemble averaged filtered sea-level 
for three-month periods at Sattahip (.56 degrees of freedom) 



spectra 
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Figure 7. Seven-year ensemhle averaged filtered sea-level spectra 
for three-month periods at Ko-Lak (56 degrees of freedom). 
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seasonal periods each contain an even number of months, with the southeast 
season ( February-March) included in the northeast monsoon season. 

The spectra from the two stations are similar in shape, as shown in 
Figure 9, for each season. Both monsoon seasons show maxima at frequency 
0.00195 cycles per hour (period 512 hours) and decrease rapidly to fre- 
quency 0.0058 cycles per hour (period 171 hours), then decrease slowly 
to frequency 0.014 cycles per hour (period 73 hours) where the energy 
tends to zero. These spectra differ in that the energy in the northeast 
monsoon season is clearly higher than that during the southwest monsoon. 

This suggests that the northeast monsoon wind has greater variation than 
the southwest monsoon. Obviously, this is the season in which sea level 
prediction most needs consideration of forcing by wind. A question which 
arises later is whether the sea level response is due to local forcing 
over the Gulf of Thailand or to remote forcing, perhaps in the South China 
Sea . 

D. METEOROLOGICAL EFFECTS 

To find connections between sea levels at Sattahip and Ko-Lak and 
meteorological processes, the relations between sea level and atmospheric 
pressure (Sattahip), sea level and the mean geostrophic wind (Bangkok- 
Saigon and Bangkok-Si ngapore) , and sea level and local surface wind (Bangkok) 
were examined. The spectrum, cross-spectrum, coherence, and phase were 
calculated by program FTFREQ (IMSL, 1980), which covered a 6-month period 
(October-March and April -September) ; the maximum number of lags was 30 
and spectral estimates have 12 degrees of freedom. 



36 



2CC 



160 



120 




Figure 8. Seven-year ensemble averaged nontidal sea-level spectra 
over the whole year (224 degrees of freedom). 
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Figure 9. Seven-year ensemble averaged filtered sea-level 
spectra according to season (112 degrees of freedom). 
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The term, "meteorological effect", generally includes the effects of 
(1) atmospheric pressure, (2) the transient wind that accompanies rapidly 
moving storms, and (3) the steadier winds on sea level (Miller, 1956). 

In this analysis, the sea level, the atmospheric pressure, and the 
wind are represented by daily values at 06 hours in the year 1962. The 
sea levels with the astronomical tide removed by low pass filter (24-24- 
25 hr.) are shown in Figure 10. The atmospheric pressures from which 
geostrophic winds were computed are presented in Figures 11 and 12. 

1 . Response of Sea Level to Atmospheric Pressure 

The change of sea-level in response to atmospheric pressure varia- 
tion in time equals the height of a column of water having the same weight 
per unit area as the change in atmospheric pressure. A pressure increase 
of 1 mb. would decrease the sea surface elevation by approximately 1 cm. 
under quasi-static conditions (Roden, 1966). 

Since the filtered sea levels and local atmospheric pressures are 
available only at Sattahip, and because of many missing values in the 
atmospheric pressure data, the response of the filtered sea levels to the 
variable atmospheric pressure were calculated from their statistics. 

The ratio of the response can be estimated formally from the ratio 

Resp = Ep/Ez 

where Ep and Ez are the variances of the atmospheric pressure fluctuations 
and the filtered sea levels, respectively. A value of 1.0 would indicate 
total response. The result for this response at Sattahip is 0.0003. 

This suggests that the atmospheric pressure does not influence sea level 
to an important extent. 
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2 . Sea Level and the Geostrophic Wind 

Geostrophic wind from daily weather records was used to infer the 
strength and direction of wind stress on the sea surface. The procedure 
is not completely satisfactory, since the relationship between surface 
and geostrophic wind is not always consistent; but within reasonable 
limits, it is useful to assume that surface wind varies linearly with 
the geostrophic wind (Petterson, 1940). 

The geostrophic winds were investigated by use of surface pressures 
at three stations (Bangkok, Saigon, and Singapore). The geostrophic 
winds found from these pressures were related to the sea levels at Sattahip 
and Ko-Lak. For example, the Bangkok-Saigon pressure difference (Figure 
1) can be viewed as representing the mean cross-gulf geostrophic wind, 
a positive pressure difference implying wind from the east. 

The spectra, cross-spectra, coherences, and phases between the . Bang- 
kok-Saigon pressure differences and the filtered sea levels at Sattahip 
and Ko-Lak are shown in Figures 13 and 14. The relationships for each 
monsoon season (October-March and April -September) at the two tide stations 
are similar. For the northeast monsoon season (October-March) , high 
coherence occurs at frequency 0.117 cycles per day (period 8.6 days), 
and the phase difference is about 0-10 degrees indicating that the geo- 
strophic wind leads sea level. For the southwest monsoon, high coherence 
occurs at frequency 0.083 cycles per day (period 12.1 days), with phase 
about 40 degrees. This seems to indicate that the northeast anomaly 
of winds produces a rise in sea level and the southwest anomaly leads 
to low sea level in the northern end of the gulf, as can be expected from 
Ekman transport. 



40 



For the filtered sea levels at Ko-Lak or Sattahip and the atmospheric 
pressure difference between Bangkok and Singapore, shown in Figures 15 
and 16, there is high coherence at frequency 0.117 cycles per day, with 
a phase difference of about 220 degrees at Ko-Lak and 0 degrees at Sattahip, 
during the northeast monsoon. For the southwest monsoon, the high 
coherence occurs at frequency 0.083 cycles per day, phase about 10 
degrees at Ko-Lak and 50 degrees at Sattahip. This shows a less consistent 
relation between sea level and the mean geostrophic wind between Bangkok 
and Singapore than for the Bangkok-Saigon mean. The implication is that 
Bangkok-Saigon pressure differences are better representative of the 
cross-gul f wind. 

3 . Sea Level Slope and the Local Mind 

The wind is one of the primary agencies that drives ocean currents 
and influences the redistribution of density within the ocean. This also 
is one of the major factors in producing sea level variations. 

The effect of local wind is most pronounced according to Groves (1956) 
as follows: 

"1) Where (or when) there are strong winds 

2) On the edge of continents; this follows from the fact that 
divergence of the horizontal wind-stress is one of the most effective 
agencies in bringing about surface slopes. A coast line will act as a 
line of divergence if the wind is blowing across it. 

3) At stations where there is an extensive shoal region adjacent 
to the tide gauge, which follows from the fact that the total horizontal 
force on a vertical column of water is independent of the water depth in 
the case of surface wind stress, although it is proportional to the depth 

in the case of the pressure gradient associated with a given surface slope". 
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The slope (by) across the northern gulf may be measured by the difference 
in sea level between Ko-Lak and Sattahip The observed local 

wind at Bangkok was examined for its relation to the slope (&>{) . The 
wind was analyzed into north-south (V^ s ), east-west (V EW ), components 
and also X(Vx), Y (Vy) components, where X is the axis along the line con- 
necting Ko-Lak and Sattahip, as shown in Figure 1. Time series of the 
wind components are shown in Figures 17 and 18. 

The spectra, cross-spectra, coherences and phases between and 
each of the four components of the local wind have been calculated and 
are shown in Figures 19-22 for each season. The cross-spectra for 
(^Y,V, vs ) and (Avj, Vy) show high coherence at frequency 0.2 cycles per 
day and phase difference about 20 degrees (Vms and Vy leadA>|) during 
the northeast monsoon. They also show high coherence at frequencies 
0.1-0.133 cycles per day, phase about 10 degrees, during the southwest 
monsoon. This means that at those frequencies the water of the upper 
gulf tends to slope downward to the east when the anomalous wind is from 
the south and vice versa, which is contrary to expectation from the 
Ekman relationship. The spectra for (A>J, Vew) and (£>J, Vx) 
show high coherence at frequency 0.117 cycles per day, phase about 
200 degrees (Vew and Vy lead^), during the northeast monsoon. They 
also show high coherence at frequency 0.083 cycles per day, phase about 
180 degrees. This is consistent with wind set-up in response to east- 
west and X components of the local wind, although this effect is secondary 
to the Ekman effect over the gulf as a whole as was concluded in Section 
V B. 
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It may be concluded that the most important relation between local 
wind and response of cross-gulf sea level slope between Sattahip and Ko- 
Lak is that in which the wind blows across the gulf and the response is 
in frequency band 0.083-0.117 cycles per day (period 8-12 days), is increased 
sea surface slope upward in the direction toward which the wind blows. 
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Figure 10. Daily filtered sea level at 
Ko-Lak and Sattahip (1962) 




45 



Figure 11. Daily sea surface pressure at 
Singapore and Bangkok (1962) 
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Figure 13. The spectrum, cross-spectrum, coherence, and phase for 
filtered sea level (Sattahip) and pressure difference (Bangkok-Saigon) , 
(12 degrees of freedom). 
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Figure 14. The spectrum, cross-spectrum, coherence, and phase for 
filtered sea level (Sattahip) and pressure difference (Bangkok-Saigon) , 
(12 degrees of freedom). 
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Figure 15. The spectrum, cross-spectrum, coherence, dnd phase for 
filtered sea level (Sattahip) and pressure difference (Bangkok-Sinqapore) 
(12 degrees of freedom). 
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Figure 16. The spectrum, cross-spectrum, coherence, and phase for 
filtered sea level (Ko-Lak) and pressure difference (Bangkok-Si ngapore) , 
(12 degrees of freedom). 
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Figure 17. Daily wind components (N-S, and E-W), wind speed 
in knots. 
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Figure 18. Daily wind components (.X and Y). , (See Figure 1), 
wind speed in knots. 
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Figure 19. The spectrum 
sea-level slope (Ko-Lak - 
(.12 degrees of freedom).. 



cross-spectrum, coherence, and phase of 
’sattahip). :and north-south wind component 
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Figure 20. The spectrum, cross-spectrum, coherence, and phase of 
sea-level slope (Ko-Lak - Sattahip) and east - west wind component 
(12 degrees of freedom). 
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Figure 21 . The spectrum, cross-spectrum, coherence, and phase of 
sea-level slope (Ko-Lak - Sattahip). and Y wind component Q2 degrees 
of freedom) . 
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Figure 22. The spectrum, cross-spectrum, coherence, and phase of sea 
level slope (Ko-Lak - Sattahip) and X;wind component (.12 degrees of 
freedom) . 
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VI. SUMMARY AND CONCLUSIONS 



Seven years of hourly sea level data from two tide stations in the northern 
Gulf of Thailand have been treated by low pass filter (using averages of 
24, 24, and 25 hours) to remove the prominent semi-diurnal and diurnal tidal 
components. Further smoothing minimized the prominent fortnightly cycle 
of spring/neap tides. It was found that the seasonal sea level variation is 
of the same magnitude as the semi-diurnal and diurnal tides (range about 5 
decimeters), and that the larger short-term sea level variations superimposed 
range up to 4 decimeters with durations of about 8 days. It was determined 
that these sea level variations are due principally to meteorological effects. 

The annual variations of sea level were studied by comparing then with 
the annual variations of atmospheric pressure and wind. The atmospheric 
pressure was found to vary directly with the sea level, and accordingly can 
account for very little of the annual sea level variation. For the comparison 
with the annual wind, sea level varies seasonally as if controlled by the Ekman 
transport which would be expected from the monsoon wind regime over the Gulf. 

The spectra, cross-spectra, coherences, and phases for the hourly sea 
levels show significant peaks at frequencies of the diurnal and semi-diurnal 
astronomical tide forces, the diurnal component dominating. The spectra 
were made free of these tidal components by low-pass filtering (24-24-25 hours) 
of the hourly tide heights. These spectra were, in turn, ensemble averaged 
over the seven-year period (Figure 8). The spectra are seen to be very similar 
for both tide stations. This implies that residual sea levels at both 
stations change in response to the same driving mechanisms. 

For the purpose of examining seasonal variation, the spectra of the 
filtered heights were ensemble averaged for each monsoon season (Figure 9). 
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The spectra show higher energy in the northeast monsoon than during the 
southwest monsoon. This is interpreted to mean that the northeast monsoon 
winds have greater variation than the southwest monsoon winds. 

The spectra, cross-spectra, coherences, and phases between the low-passed 
sea level and the geostrophic wind were also computed (Figures 13-16). 

These show that the sea levels are in phase with the geostrophic winds computed 
from pressure measurements at Bangkok and Saigon. This suggests that sea 
levels relate to fluctuations in the monsoon winds in a way consistent with 
piling up of water at the northern boundary of the gulf by northward Ekman 
transport and vice versa for southward transport. 

The relationship of the low-pass filtered sea level slopes between the 
two tide stations, located on opposite sides of the northern gulf, and the local 
wind was also studied. The cross-spectra (Figures 19-22) show high coherence 
in the frequency band 0.083-0.117 cycles per day for the east-west and X- 
components of the wind. This indicates that in this frequency band the sea level 
slopes upward across the gulf in the same direction as that toward which the 
anomalous wind blows, the response being most coherent for periods of 8-12 days. 
This effect is clearly secondary to the effect of Ekman transport in producing the 
seasonal sea level variations observed in the northern Gulf of Thailand. 

It is clear that the sea level variations studied here are affected by 
the monsoon wind and its variations. The observed sea level variations are not 

totally explained by wind-forcing over the gulf. They are possibly, in part, 
responses to forcing over the South China Sea. Future studies should 
examine: (1) correlation between gulf sea levels and wind forcing in the 

South China Sea and (2) dynamical /numerical models for responses to local 
forcing and to forcing over and from the South China Sea. 
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APPENDIX A 



RAW AND FILTERED HOURLY SEA LEVEL DATA 

The following time-series graphs for Sattahip and Ko-Lak show: 

(1) raw or unfiltered hourly sea levels (lower), (2) semi-diurnal and 
diurnal hourly water levels extracted using a 24-24-25 hour filter (upper), 
and (3) the residual hourly water levels remaining (middle). The water 
levels are in decimeters and measured from OOOOZ on 1 January. Each 
graph covers one year of data. 
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APPENDIX B 



HOURLY SEA LEVEL DATA FILTERED 
WITH GODIN AND 28-DAYS FILTER 

The following time-series graphs for Sattahip and Ko-Lak 
show: (1) residual hourly water levels remaining after application of 

Godin (24-24-25 hour) filter to the raw hourly data (upper graph; also 
shown in Appendix A), and (2) the same data after further smoothing 
using a 672-hour (28-day) filter. Water levels are in decimeters and 
the time scale is in hours measured from 0000Z on 1 January. Each 
graph covers one year of data. 
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APPENDIX C 



SPECTRA, CROSS-SPECTRA, COHERENCE, 

AND PHASE OF RAW SEA LEVEL DATA 

The following graphs for Sattahip and Ko-Lak show (from the 
top): (1) phase, (2) coherence, (3) cross-spectra, and 

(4) spectra for three-month period (8 degrees of freedom). 
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APPENDIX D 



SPECTRA, CROSS-SPECTRA, COHERENCE, AND 
PHASE OF LOW-PASS FILTERED SEA LEVEL 



The following graphs for Sattahip and Ko-Lak show 
(1) phase, (2) coherence, (3) cross-spectra, and 
three-month periods (8 degrees of freedom). 



(from the top): 
(4) spectra for 
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